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y—{ Abstract. The optimisation of two-dimensional (2D) lattice ion trap geometries for 

trapped ion quantum simulation is investigated. The geometry is optimised for the 

(— ! highest ratio of ion-ion interaction rate to decoherence rate. To calculate the electric 

P^ I field of such array geometries a numerical simulation based on a "Biot-Savart like 

law" method is used. In this article we will focus on square, hexagonal and centre 

rectangular lattices for optimisation. A method for maximising the homogeneity of 

^ trapping site properties over an array is presented for arrays of a range of sizes. We 

(sJ show how both the polygon radii and separations scale to optimise the ratio between 

the interaction and decoherence rate. The optimal polygon radius and separation for a 

CNJ 2D lattice is found to be a function of the ratio between rf voltage and drive frequency 

^' applied to the array. We then provide a case study for ^^^Yb"*" ions to show how a 

r two-dimensional quantum simulator array could be designed. 

(N 

en 

O 1. Introduction 

L| Trapped ions possess long lived addressable internal states and can be highly decoupled 

*0 from their environment. This makes them an important tool in the development 

;_( of quantum information processing [H [2] and quantum simulation [31 IH El El E]- 

When used for quantum simulation they enable complex spin systems, among others, 
to be investigated beyond the practical limitations of classical computation. For 
example, trapped ions have been used for quantum simulations of the evolution of 
paramagnetic into (anti-)ferromagnetic order in a spin system, [HI [9] and frustrated 
anti-ferromagnetic Ising interactions, [TUIIII]. These first simulations were carried out 
using one-dimensional trapping arrays and state dependant forces applied using laser 
beams. 

More complex simulations will require ions trapped in 2D arrays and interaction 
schemes compatible with these. Advances towards 2D trapping arrays suitable for 
quantum simulations have been made by trapping ions in a millimetre-scale mechanically 
fabricated metal mesh [12] and by the successful implementation of microfabrication 
techniques for ion traps [13]. In addition, interaction schemes based on large oscillating 
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or static magnetic field gradients have been proposed [HI [15] and demonstrated in 
an on-chip microwave gate p^ [17]. With the recent advances in the field, ion trap 
quantum simulations using large scale 2D ion trap lattices are within the reach of current 
technology. 

In order to create a 2D array of trapped ions a repeating 2D surface geometry 
is required. Decoherance due to anomalous heating is a major issue for large scale 
quantum simulations. As this heating scales approximately as r~^ [TH], where r is the 
ion height above the trapping surface, it is advantageous for ions to be trapped high 
above the surface. However, when individual surface microtraps are placed together 
so that their separation is less than around twice the ion height the individual electric 
fields start to overlap and distort the resulting trapping fields [19] . In extreme cases this 
can lead to the traps combining to produce a singular trapping zone. To compensate 
for this the electrode structure has to be altered when operating within this regime 
^19j . Schmied et al [19] have investigated surface-electrode geometries and developed 
an algorithm that optimises geometries to maximise the electric field curvatures of 
individual trapping sites for arbitrary ion heights and separations. Individual trapping 
sites shown in [T^ [20] were optimised using this algorithm leading to non-intuitive 
electrode patterns which can contain many isolated radio-frequency (rf) and static 
voltage electrodes. Another proposal [21] working outside this regime uses rf electrodes 
with controllable rf voltages to lower trap frequencies and decrease ion-ion distances 
and, therefore, increase interaction strengths. However, this requires the use of multiple 
independent rf electrodes and individually controllable rf voltages posing an additional 
experimental complication. 

In this paper we present an optimisation process for ion trap topologies based on a 
single island of rf electrode, reducing the requirement for buried rf wires and multiple 
rf electrodes. We focus on the development of an optimum lattice geometry where the 
ratio of coupling rate to the decoherence rate due to ion heating is maximised and made 
homogeneous across the lattice. This is achieved by minimising the secular frequencies of 
the trapping sites whilst, simultaneously, keeping the trapping depths above a minimum 
trap depth (for illustrative purposes we use O.leV) to allow for successful operation of 
the proposed 2D lattice designs. We will concentrate on the optimum lattice topology 
for hexagonal, square and centred rectangle lattices. An investigation is also carried out 
on how optimal geometries depend on the overall lattice size, and we discuss the choice 
of and scaling for experimental parameters such as rf voltage, drive frequency, ion mass 
and electric field noise density. 

Additionally it is possible that two-dimensional ion arrays of this type could also 
be used for quantum computation. For example, cluster state quantum computation 
could be carried out in such a system [221 [23] • However, additional constraints may also 
have to be taken into account. 
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2. Ion-ion interactions and Lattice Geometry 

2.1. Ion-ion interactions 

Two-dimensional lattices of ions can be used as a quantum simulator for many body 
spin-1/2 systems [H]. Forces such as the trapping potential, Ft = —muf, and 
the coulomb force, Fc = — e^/(47reoA^), between the individual ions determine the 
equilibrium position. Laser beams or magnetic field gradients can be used to impart 
an additional force to the ions which displace the ion(s) depending on its internal 
state. This displacement leads to a change in the coulomb force and thereby displaces 
the neighbouring ion(s) dependent on their own internal state. This state dependant 
coupling is given by 



3F^ 

J=^ (1) 

moo'' 

where m is the mass of an ion, F is the magnitude of the state dependant force applied 
to each ion and u is the trap's secular frequency. We will consider how this force is 
applied later in the article. /3 is the ratio of the change in the Coulomb force to the 
change in the restoring force due to the displacement of the ions caused by the state 
dependant force and is given by [23] 

^ " 27reomw2^3' (2) 

where A is the ion-ion spacing. There are two cases to consider, when /5 > 1 the change 
in the Coulomb force, SFc, due to the displacement of an ion is dominant over the 
change in the restoring force, 6 Ft- This results in an interaction over a large number of 
trapping sites. When /3 < 1, the opposite is true resulting in an interaction which decays 
rapidly across the array. Trapping ions in a 2D array of microtraps makes it possible 
to satisfy the condition that (3 < 1 allowing systems with short range interactions to be 
simulated. An illustration of this system is shown in figure [TJ 

It is important to consider sources of decoherence when designing a 2D ion trap 
array. The internal state of an ion can remain coherent for lO's of seconds [251 126] . 
However, motional decoherence due to anomalous heating of ions will be an important 
factor during quantum operations within small scale ion traps as the implementation of 
spin dependant couplings involves the use of motional states of the ion. The coupling, J, 
will be observable if the coupling time, Tj = 1/ J, is less than the motional decoherence 
time in the system and, therefore, the ratio of these two times is an important parameter 
of the system and is given by 

Tn 



where HH 



AmuH , 



Optimisation of two-dimensional ion trap arrays for quantum simulation 
(a) 




(b) 




Figure 1. Pictorial diagram of three ions in adjacent traps. The ions feel a Coulomb 
force indicated by the springs between each ion and can be placed in two different 
states indicated by their colour, (a) Pictorial diagram showing the case with no state 
dependent force present, (b) Pictorial diagram showing how the system reacts to the 
presence of a state dependent force, F. Here the ions feel a change in the Coulomb 
force, 6Fc, due to the displacement of the ions and a change in the restoring force, 
SFt. 



Here S'£;(a;) is the electric field noise density [ISl [IS]- In order for an interaction to occur 
on faster time-scales than the decoherence in the system, we require Kgim > 1 and it is 
the aim of the optimisation process presented in this work to optimise the geometry in 
order to maximise this parameter. To acquire an understanding of how a geometry can 
affect Ksim it is necessary to determine its form with respect to the geometry variables. 
The form of Tj can be found by substituting equation [2] into [T] and is given by 



Tj 



e2^2 



The Kgijn parameter can then be expressed as 



(5) 



(6) 



S'£'(cij)7reomC(;'^y43 
The secular frequency, u, of a trapped ion [27] can be expressed as a function of a 



defined as 



a 



V 



(7) 



where V is the amplitude of the rf voltage applied to the trap and Q is 27r times the 
drive frequency in Hz, yielding. 



U 



y2mVLr'^ y2mr'^ 



(8) 



where r is the height of an ion above the surface, e is the charge of an electron and r]geo 
is an efficiency factor which can range between zero and one depending on the form of 
the geometry 
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The secular frequency given in equation M can then be used along with Se{(^ 



'=:»,-4 



r~ u , where S is a coefficient that can be experimentally obtained and depends on 
the temperature and surface of the trap electrodes (see fT3] for a listing) to re-express 
equation [6] in the form 



K. 



Sa^^aeo^reo^^ ' 



(9) 



To further understand how the geometry effects the K^jm value we will now introduce 
the parameters of a lattice geometry and relate them to equation |9j 

2.2. Two-dimensional ion trap lattice geometry 

A lattice is a regular tiling of a space by a primitive unit cell. Previous works [I2l |2T] 
concentrate on lattices created from square unit cells. In total there exist five types of 
cell which can be used to form a 2D lattice: centred rectangular, hexagonal and square 
as shown in figure [2| and rectangular and oblique [28] . The rectangular and oblique 
structures are not considered in this work due to their non-uniform ion-ion distances. 

Figure |2] shows the polygon-polygon separation which is equal to the ion-ion 
distance, A, in equation |9} The polygon radius, R, along with the separation, will 
determine the height above the surface at which the ion is trapped, r, with larger 
polygon radii yielding higher ion heights. Another variable to be considered is the gap 
between the outer polygon in the array to the edge of the rf electrode, g. This can 
be used to alter the homogeneity of the individual trapping sites within the array. In 
general a non-homogeneous system results in spin dependant coupling rates which are 
a function of the lattice site, posing a significant problem for the scalability of such an 
array [29]. 
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A 








R 
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Figure 2. Diagrams showing the polygon radii R, the separation between the polygon 
centrers, A, and the distance between the last polygon and the edge of the rf electrode 
(shown in grey), g. (a) Diagram showing a three by four ion trap surface array 
consisting of six sided polygons arranged with square unit cells, (b) Diagram showing 
a similar surface array arranged into hexagonal unit cells, (c) Diagram showing a 
surface array arranged into centred rectangular unit cells. The unit cells are indicated 
by dashed lines. 
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3. Simulation of lattices 

To determine the electric field produced by a two-dimensional array a method based on 
the Biot-Savart like law described by Oliveira and Miranda [30] was used. This method 
calculates the electric field produced by an arbitrarily shaped two-dimensional electrode 
which is held at a potential V whilst the rest of the plane is held at a potential of zero. 
The electric field observed at a given point, X, in space due to such an area held at a 
potential and bounded by a path C is given by [30] 



^ V f (x — x') X ds , , 

^ ' zn J(j \x-x\ 

where the curve, C, bounds the electrode and x' and x are vectors that locate the source 
point and field point respectively. By calculating the electric field in this manner an 
assumption is made that there is no gap between the areas held at the potential, V, and 
the areas held at zero. In microfabricated surface traps, gaps between the electrodes 
are required and typically range from 3^m - lO/xm [13]. If, however, these gaps are 
small in comparison to the electrode structures they will not alter the trapping fields 
significantly [20l|3Tl[32]. The electric fields of individual electrodes can then be combined 
to determine potential nils and, therefore, trapping positions in the 2D trap arrays, using 
the numerical Gauss-Newton algorithm. The secular frequencies, trap depths and ion 
heights at these positions can then be determined. 

To calculate the error of our numerical integration we compared simulations of five 
wire symmetric surface trap geometries with different central static electrode and rf rail 
widths in the gapless plane approximation using the method of Oliveira and Miranda 
[30j to results obtained with analytical equations described by House [31]. In all the 
geometries simulated the two rf rails were of equal width. Similarly to House, the outer 
static voltage electrodes were approximated as an infinitely long ground plane although 
the length of the inner rails were set to 3000 /im instead of infinite. A selection of these 
simulation results are shown in table [U 

In these results a general error for the ion heights and secular frequencies of less 
than 2% and 3% respectively was found, which leads to a maximum error in Kgi^ of 
10%. For the following simulations it is therefore assumed that the maximum Kgim error 
is 10%. Additionally, numerical simulations of the geometries were carried out using 
methods described in [53], which indicate similar errors and trends for the ion height 
and secular frequency as the results obtained with the Biot-Savart like method. 

4. Lattice geometry optimisation 

In this section we show how the parameters of the lattice geometry (discussed in section 



2.2) can be optimised to achieve the highest possible K^jm value across the array for a 
given set of experimental parameters. To do this, we first show how to maximise the 
homogeneity of individual site properties over an array by varying the distance between 
the outer polygon in the array to the edge of the rf electrode, g, and show how this scales 
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Electrode parameters 


Simulations 


House 


equations 


rf width 

[/im] 


central static 
electrode width 

[/im] 


rf Volt. 

[V] 


rf freq. 

[MHz] 


r 

[/im] 


[MHz] 


r 

[/im] 


u 
[MHz] 


100 


50 


250 


75 


55.8 


6.86 


55.9 


6.87 


100 


50 


500 


60 


55.8 


4.29 


55.9 


4.29 


200 


100 


250 


30 


110.1 


2.17 


111.8 


2.20 


200 


100 


500 


40 


110.1 


3.26 


111.8 


3.30 


500 


150 


250 


20 


165.4 


1.48 


167.7 


1.52 


500 


150 


500 


25 


165.4 


2.37 


167.7 


2.43 



Table 1. Table showing the secular frequency, w, and ion height, r, for different five 
wire surface trap geometries as calculated by the analytical method in House [31 and 
simulated by the method used in this work based on the Biot-Savart like law [3D]. 



with lattice size. We will then use the homogeneous arrays to calculate the optimum 
number of sides, n, a polygon within the array should possess in order to maximise Kgim- 
We then outline a method for the optimisation of the polygon radii, R, and separation, 
A, of an array and show how these vary with increased lattice size and ion mass. 



4.I. Increasing the homogeneity of Kgi^ across the array 

This is achieved by ensuring homogeneous secular frequencies, ion heights and trap 
depths across all the array sites. As shown in figure [3] the K^jm of trapping sites in an 
array can be altered to approach a common value if the distance, g, between the edge 
of the outer polygon and the edge of the rf electrode containing the polygon array is 
adjusted. As the value of g is increased, the Kgim value of the sites towards the centre 
drop, however, the outer sites Ksjm value rises, resulting in the properties converging 
towards each other. If the distance g is increased further beyond the point at which 
maximum homogeneity occurs the outer site properties drift away from those of the 
central sites and, therefore, decrease the homogeneity of the array. 

To provide a value of g which is universal for all lattice sites, its value is given in 
units of lattice side length, L. The lattice side length is determined by the polygon 
separation, A, and radius, i?, and can be expressed as 



L = {M -l)A + 2R, 



'111 



where M is the number of lattice sites along one side of an array. 

In order to quantify the arrays homogeneity, H is defined as the average deviation of 

of the central site and is given by 



Ksjm of each lattice site from the K 



H 



1 ^ 

-Y 

n=l 



1- 



K ■ 



K, 



(12) 



where iV is the total number of trapping sites in the lattice. 
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Figure 3. Diagram sliowing tlie effect of varying tlie distance g on the scaled Kgim 
value of the individual trapping sites. The Kgim values shown are scaled with that of 
the central site, (a) Representation of a 5 by 5 square type lattice array indicating the 
axis labelling, (b) Slices across the array (indicated by the dotted line in (a)) for g/L 
values of 0.1, 0.2, 0.5 and 1. 



Lattice Type 


a 


b 


B 


Square 


0.20±0.01 


5.21±0.38 


0.74±0.03 


Hexagonal 


0.39±0.08 


3.76±1.40 


0.54±0.14 


Centre Rectangular 


0.31±0.04 


7.84±2.43 


0.76±0.11 



Table 2. Table showing a and b values for the fits which describe g/L as a function 
of the number of sites in the lattice. 



Figured shows H for a five by five square type unit cell lattice for < g/L < 1.5. 
The maximum homogeneity, and thus the optimum g/L, is found when H is minimised. 
The error associated with H is given by 



(^H 



\/N{aKsim? + N{aKs 

N 



(13) 



where the error on all ¥^sim values is 10%, as shown in section [3] This yields an overall 
percentage error on H of 0.13/viV%- 

Figure p^ shows the optimum g/L for hexagonal, central rectangular and square unit 
cell lattices of different sizes. The curves are found to be described by an equation of 
the form g/L = a + bN~^ with a, b and B values for different lattice types shown in 
table [2] For large lattices, g/L is independent of N, as trapping sites close to the edge 
of the lattice are influenced only by the electric field created by that edge. In small 
lattices, however, the optimum g/L increases as the effect of the electric field from the 
opposite edge of the lattice increases. 
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g/L 

Figure 4. Graph showing the average deviation of the Kgim of each lattice site from 
the Kgim of the central site, H, for a five by five square type unit cell lattice for 
< g/L < 1.5. The error on H is given by ^^j^H and the error of the minimum of 
H is determined by observing the spread of g/L which agrees, within error, with the 
minimum position. 
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Figure 5. Graph showing the optimum g/L as a function of the total number of 
trapping sites, N, for square lattices (square markers), hexagonal lattices (circular 
markers) and centre rectangular lattices (diamond markers). The curves are given by 
g/L = a + bN~^ with a, b and B values for different lattice types shown in table |2] 
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4-2. Optimising the number of polygon sides 

We now investigate the optimum number of polygon sides, n, providing the highest Kgim 
value on the central trapping site (located above the central polygon) of a lattice. 

To ensure the results are universal for all lattice geometries, g/L is set to the value 
which maximises the homogeneity of each array, and all other parameters are scaled by 
normalising Kgim to that of an identical geometry with polygons of 100 sides. This also 
allows comparison between the different types of lattices. 




20 30 

Number of polygon sides 



50 



Figure 6. Graph showing the relationship between the number of polygon sides and 
Ksim for square (square markers), hexagonal (circular markers) and central rectangular 
(diamond markers) unit cell lattices. The dashed lines show the asymptotes of 1 and 
0.95 of the scaled Ksim value. 

Figure |6] shows the scaled Kgim for the central site as the number of polygon sides is 
varied. As the number of sides is increased the value of K^jm approaches an asymptote, 
shown by the upper dashed line. This indicates that the best geometry will be made 
from circular electrodes. However, simulation times grow with increasing polygon side 
number and so it is advantageous to reduce this number to a minimum. It is shown 
that ~95% of the maximum achievable Ksim (indicated by the lower dashed line) can 
be achieved with around >25-30 sides in the polygons. 



4-3. Optimisation method for polygon separation and radius 

In this section we now maximise the Kgim of any arbitrary geometry. We will then go 
on and determine optimum geometries and show how they scale and, ultimately, are 
determined by a = V/Q. We use the value of g determined in section 4.1 in order to 
provide the maximum Kgim homogeneity across the lattice, and set the number of sides 
to 25 as this provides a good approximation to the optimum circular geometry while 



keeping the simulation time at a minimum, as shown in section 4.2 
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When considering any fixed arbitrary geometry, equation |9] shows that Kgim can be 
maximised by reducing the value of a. However, the minimum achievable a is hmited 
by the lowest usable trap depth, as the trap depth is proportional to a^ and is given by 

Tn = ^-^ (14) 

where e is the charge of an ion and C is a geometrical factor which is a function of A 
audi? [3I]. 

We will now focus on finding optimal geometries which we define as geometries, 
which yield the highest value of Ksim for a given value of a. This will be carried out 
by fixing the trap depth at a reasonable minimum value (we use 0.1 eV for illustration 
purposes) which, as discussed above, provides the maximum K^jm for a given geometry, 
and then investigating the dependence of polygon radius, separation and ion height 
with a. To determine these dependencies a K^j^ contour plot is made by calculating 
the Ksim over a range of polygon separation. A, and radius, R, with a resolution of l/im. 
The range of polygon radius and separation used should not create traps with inter-well 
barriers of less than the minimum trap depth value, and to ensure this the polygon 
radius was kept to less than a third of the polygon separation. For each combination of 
polygon separation and radius a value of a is found which yields the minimum trap depth 
and, thus, maximises the Ksim of the particular geometry. By following this method 
one can obtain the a required to achieve the minimum trap depth, the ion height, r, 
and Ksim for each geometry. From the resulting data the polygon separation and radius 
which yields the highest Ksim for a given a (the optimum geometry) can then be found. 
A graphical example of such data is shown in figure [7j 

It can be seen from this method and the example data in figure [7] that the highest 
Ksim will be achieved with an infinite value of a, R and A. However, other effects may 
limit the magnitude of a. In order to determine a limit on a it is, therefore, necessary 
to describe the various array and trapping field dependant properties (such as secular 
frequency, ion height and Ksim) in terms of a. 

By plotting these optimum parameters (polygon radius, separation and ion height) 
as a function of a, as shown in figures [sta), (b) and (c), linear relationships of the form 



and 



r = kra (15) 

A = kAa (16) 

R = kna (17) 



are found for the optimal geometries. The values of kr, kA and k^ are dependant on the 
number of trapping sites in an array, as shown in figures [9ta),(b) and (c) respectively, 
for lattices made from square type unit cells of polygons. 



It is important to stress that equations 15, 16 and 17 are only valid in the case of 



optimal geometries, which depend solely on a. With this in mind, it is now possible to 
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Figure 7. Example graph sliowing liow tlie Kgim (absolute values indicated by 
numbers on contour lines) varies as a function of polygon radius and separation. The 
graph also indicates that the value of a increases as the radius and separation are 



increased. This data was obtained using the method described in section 4.3 with a 
polygon separation and radius resolution of 1/im and a minimum trap depth of 0.1 eV. 
The value of a in the figure ranges from zero to sal. 5 VMHz~^. The impossible region 
describes geometries where individual trapping sites start to combine to a single one 
and so posses a polygon radius, R, greater than or equal to a third of the polygon 
separation, A. 



100 




o;[VMHz] 
(c) 
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a [V MHz] 




0.5 1.0 

or [V MHz] 



Figure 8. Graphs showing the ion height (a), polygon separation (b) and polygon 
radius (c) of an optimised lattice as a function of the ratio a. In all cases the plots are 
shown using a = V/fl where fl is 2tt times the drive frequency in Hz, and for arrays 
made from square type unit cells of polygons with 81 sites and for ^^^Yb+ ions. 
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Figure 9. Graphs showing the value of kr (a), kA (b) and ku (c) as a function of tlie 
number of trapping sites, N. In all cases the plots are shown using a = V/fl where fl 
is 27r times the drive frequency in Hz and are for arrays made from square type unit 
cells of polygons using ^^^Yb+ ions. 



re-express the secular frequency in equation [8] to describe the secular frequency of an 
optimised geometry as: 



U 



V2k; 



ma 



Ksim in equation [9] can also be re-expressed to describe that of an optimised 



geometry: 



K. 






(19) 



Equation 19 shows that, for optimal geometries, Ksim is proportional to a (as the 



value of a determines the electrode dimensions to be used) and so, to produce an array 
with a high Ksim for a given number of lattice sites (as kr and kA are a function of 



the number of trapping sites), a large value of a is preferable. Equations 16 and 17 
show that the optimum geometry size is proportional to a. It, therefore, follows that 
larger lattice geometries will produce larger values of Ksim- This effect is illustrated 
in figure [7] which shows the K^j^ as a function of polygon radius and separation. For 
optimised lattices, the optimal radius and separation will fall on a line described by 
^ = (^a/^k) R, with higher values of a required for higher values of separation and 



radius as shown in equations 16 and [T7| respectively. 
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The heating rate in ion traps has a strong dependency on the ion height (oc l/r~'^) 
|18j . A large K^j^ is achieved with large values of a, resulting in large ion heights, 



as shown in equation 15 It, therefore, can be concluded that a different scaling of 
the heating rate (for example in cryogenic systems) does not change the optimisation 
process and optimal geometries. 

It has now been shown that the optimal geometry for a given minimum trap depth 
and ion mass is determined solely by the value of a. 

Optimal geometries and their Ksim values (in units of 1/a^) can now be found by 
creating contour plots (such as shown in figure [T]) for different values of lattice size, N, 
lattice unit cell type and ion mass, m. The error on the Kgim value, calculated from 
equation |9| was determined to be ±10% by comparing the secular frequency and ion 
heights obtained using the program with those predicted by House's analytical solutions 
for a five wire surface trap geometry |31j . 

5. Optimisation results and analysis 

In this section, optimum polygon separations. A, and radii, R, are obtained using 
the method outlined above for square, hexagonal and centre rectangular unit cell 
type lattices. These are shown as function of lattice sizes and ion masses with the 
experimental constraint, a, scaled out. Throughout this optimisation example, ^^^Yb"*" 
ions will be used unless otherwise stated. 



Figure 10 shows how the optimum scaled radius, R/a, and separation. A/a, of 
polygons vary as a function of lattice size for ^^^Yb"*" ions. As explained in the previous 
section we have assumed a minimum trap depth of 0.1 eV for illustrative purposes. 
It can be seen from this figure that as the size of the lattice increases, the optimum 
polygon radius and separation asymptotically tend towards values representative of an 
infinitely large lattice. This is expected as once a lattice becomes large enough, the 
addition of extra lattice sites will represent only a small change in the overall electrode 
geometry and, therefore, produce a small change in the electric field produced by the 
geometry. When the lattice is small however, additional lattice sites will represent a 
larger change in the geometry and will, therefore, cause a larger change in the electric 
field. Figure 11 shows how the scaled Ksim/{F'^a^) scales as a function of the number 



of lattice sites, N, using scaled optimum polygon radii, R/a, and separations. A/a. 
The state dependant force F will be considered in more detail in section |6.2[ Due to 
the dependency of Ksim on the geometry, the relationship between Kgim/iF'^a^) and the 
number of sites is expected to be of similar form to that for optimal polygon radii, R/a, 
and separation. A/ a, with the maximum KsJm/(F^tt^) asymptotically tending towards 
a value representative of an infinitely large lattice. 



Using the data shown in figures 10 a), [lO[b) and 11 the optimum radii and 



separation of the polygons were found to follow a c + dN~^ and / + gN'*^ relationship, 
respectively, while the maximum Ksim/(F^a^) follows a k + IN~^ trend. The values of 
c, d, E, f, g, G k, I, and Q are shown in tables |3} |4]and|5| 
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Lattice Type 


c 


d 


E 


Square 


-5±1 


101±7 


0.29±0.03 


Hexagonal 


1±3 


85±6 


0.40±0.08 


Centre Rectangular 


-2±2 


68±6 


0.34±0.05 



Table 3. Table showing c, d and E values for the fits which describe R/a as a function 
of the number of sites in the lattice. 



Lattice Type 


f 


g 


G 


Square 


-136±7 


457±5 


0.15±0.01 


Hexagonal 


13±24 


547±288 


0.48±0.21 


Centre Rectangular 


40±15 


831±282 


0.57±0.13 



Table 4. Table showing / and g and G values for the fits which describe A/a as a 
function of the number of sites in the lattice. 



Lattice Type 


k 


1 


Q 


Square 


-(2.69±2.97)xl034 


(3.09±L31)xl036 


(0.61±0.13) 


Hexagonal 


(0.52±6.00)xl033 


(3.51±1.54)xl036 


(0.86±0.12) 


Centre Rectangular 


-(3.09±6.12)xl035 


(5.23±5.11)xl037 


(0.65±0.26) 



Table 5. Table showing k, I and Q values for the fits which describe Ksim/(P'^ct^) as 
a function of the number of sites in the lattice. 




(b) 
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Figure 10. (a) Graph showing how the optimum polygon radius, R/a, varies as 
a function of the number of sites, (b) Graph showing how the optimum polygon 
separation, A/a, varies as a function of the number of sites. For both (a) and (b) the 
results shown are for square (square markers) , hexagonal (circular markers) and centre 
rectangular (diamond markers) unit cell lattices with ^^^Yb"*" ions. 
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Figure 11. Graph showing how the optimum Ksim/(P'^ct^) varies as a function of the 
number of sites for optimum lattices with ^^^ Yb"*" ions. This is shown for square (square 
markers) , hexagonal (circular markers) and centre rectangular (diamond markers) unit 
cell lattices. Here F is a state dependant force applied to the ions in the lattice. 



Lattice Type 



o 



P 



Square 

Hexagonal 

Centre Rectangular 



-56±6 138±10 

-34±7 88±11 
-48±3 110±6 



Table 6. Table showing o and p values for the fits which describe R/a as a function 
of ion mass. 



Using the data shown in figures [l2[a) and ^2[b) the optimum radii and separation 
of the polygons was found to follow a o-\-pm~^'^ and g+sm~°'^ relationship, respectively, 
with values of o, p, q and s shown in tables [6] and [7j Figure 13 shows how the optimum 
KsJm,/(F^a^) varies as a function of the mass of the trapped ion, m, for 220 (square type 
unit cells) and 225 (hexagonal and centre rectangular type unit cells) trapping sites. It 
is found that the optimum KsJm./(F^a^) scales as -u + vm~^-^, with the values of u and 
V shown in table [HI 



We note, as shown in figure 12 , that as the mass of the ion is increased, the polygon 



radii and separation will have to be decreased in order to provide trapping regions with 
a depth of above 0.1 eV for a given a for 220 (square type unit cells) and 225 (hexagonal 
and centre rectangular type unit cells) trapping sites. It is also clear to see that ions 
with lighter masses will provide higher Ksim/{F'^C(^) values but will require larger lattice 
geometries compared to heavier ions. 
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Lattice Type 


q 


s 


Square 


-176±47 


489±65 


Hexagonal 


-311±15 


678±22 


Centre Rectangular 


-402±28 


911±51 



Table 7. Table showing q and s values for the fits which describe A/ a as a function 
of ion mass. 




so 100 150 

Mass of trapped ion [amu] 



200 



50 100 ISO 

Mass of trapped ion [amu] 



200 



Figure 12. (a) Graph showing how the optimum polygon radius varies as a function of 
the ion mass for 220 (square type unit cells) and 225 (hexagonal and centre rectangular 
type unit cells) trapping sites, (b) Graph showing how the optimum polygon separation 
varies as a function of the ion mass. In both graphs this is shown for square (square 
markers) , hexagonal (circular markers) and centre rectangular (diamond markers) unit 
cell lattices and the polygon radii and separations are scaled with a. 



3.0x10^* 



'^ 2.5x10 



< 2.0x10^** 




100 150 

Mass of ion [amu] 



200 



Figure 13. Graph showing how the optimum Ksim/(P'^ce^) varies as a function of the 
ion mass for 220 (square type unit cells (circular markers)) and 225 (hexagonal and 
centre rectangular type unit cells (square markers and diamond markers respectively)) 
trapping sites. 
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Lattice Type 


u 


V 


Square 


-(5.25±3.00)xl035 


(7.96±1.72)xl036 


Hexagonal 


-(1.25±0.60)xl035 


(5.08±1.29)xl036 


Centre Rectangular 


-(2.56±1.28)xl034 


(3.27±0.86)xl036 



Table 8. Table showing u and v values for the fits which describe Ksim/{F'^a^) as a 
function of the ion mass. 



6. Constraints on a 

In this section we will discuss the considerations which could limit the value of a. To 
do this we will show how the power dissipation in a chip trap, the quantum simulation 
error and the interaction time vary as a function of a. This is important as from this a 
value of a can be determined for a given experiment, which will be shown in section [7j 

6.1. Power dissipation in optimised arrays 

The power dissipation of an ion trap chip is determined by the voltage, V, frequency, 
Q, as well as the capacitance and resistance of the trap itself. This may, for a given 
capacitance and resistance, affect the value of a (the ratio between the rf voltage and 
drive frequency) which can be used. As the value of a is used to determine the optimum 



polygon radii and separation of a geometry, as shown in figure 10 , for a given number of 
sites, A^, and stability parameter, q, it is important to know how the power dissipation 
varies as a function of a. 

The power dissipation of a chip is approximately given by [I3] 

1, 



Pi 



D 



-V'C'RVt' 



(20) 



where C and R are the capacitance and resistance of the chip. It is not possible to apply 
any combination of V and f2 to a geometry as they must be chosen so that the ion is 
stably trapped with a stability parameter given by [121 EZ] 



Q 



2er]geoV 2er]geoa 



(21) 



mr'^Vt'^ mr'^Vt 

between zero and 0.9, where e is the charge of an electron. 

The ion height, r, of ions trapped in the optimised lattices shown in this work have 
been found to be linearly proportional to a. This relationship is shown in figure IS] (a) 
with the constant of proportionality, kr found to be k. 60.7 mV^^s^^ for the example 
case of square type lattice with 81 sites using ^^^Yb^ ions. Considering one particular 
ion height, tq, and substituting for tq 



On 



mk'^aoq 



krao and rearranging equation |2T] for Q yields 

(22) 
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This equation can be re-expressed for Vq by noting that Vq 



fioao: 



Va 



mk'j.q 



(23) 



Equations 22 and 23 show that, for a given ion mass, m, ion height, tq, stabihty 
parameter, g, and number of trapping sites in the array (as k^ is a function of the number 
of trapping sites), there is one unique voltage, Vq, and unique parameter ao- This means 
that a given ion height (and, therefore, a chosen value of a) determines both the voltage 
and drive frequency to be applied to the trap. 



To express the power dissipation, P^, in terms of a equation 22 and 23 can be 



substituted into equation [21] giving 



D 



Im^q'^a^ 



(24) 



Equation 24 shows that as a is increased, the power dissipated is reduced. This means 



that the power dissipation is low for high values of a and, as high values of a provide high 
values of ¥^sim (see figures 11 and 13), power dissipation will not impact on producing 
high values of 'Ksim in optimised geometries. However, a low value of a will result in a 
high power dissipation in the chip and, so, the maximum allowable power dissipation in 
a chip will determine the lowest a which can be applied to a geometry. 



6.2. Quantum simulation error 

An upper limit on the value of a can be obtained from an estimation of the error of a 
quantum simulation using the method described in f3^, where the error for the Ising 
model is given by 

1 



Eo ^ -v' E(2^ + 1) i[[0{t), ^'M ' ^M > • 



(25) 



Here n is the mean radial phonon number of the ions, O is the observable of the quantum 
simulation and ?7 is a parameter which characterises phonon displacement caused by the 
state dependant force and is given by [3i 



V 



F^/m 



niuj] 



hcu 



(26) 



where m is the mass of a trapped ion and uj/2tt is its secular frequency. 

If O is an M-site observable then there exist M non-vanishing commutators (for 
example a two-site correlation function (M = 2) or a spin average (M = 1)) and so the 
error on the simulation will not be dependant on the number of ions, N, in the array 



The error in equation 25 can now be re-written as 



En 



F^M{n 



2hmuj^ 



(27) 



Equations ^9] and 27 show that both the K^j^ and the error of the simulation, Ei 



0, 



are proportional to the square of the state dependant force, F. It follows that the way 
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in which this force is apphed to the ions will determine the dependence of the simulation 
error on a. A discussion on the possible effects of the heating rate on the error can be 



found in Appendix A In this work we will consider applying this state dependant force 
via laser beams and magnetic field gradients. 

To calculate the laser power required to achieve a force, F, it is assumed, for 
illustrative purposes, that the laser beam is focused to a sheet given by 25 /xm multiplied 
by the width of the array. The laser intensity required to provide a state dependant 
force, F, can be provided by a laser beam of power, P. If the output power of the laser 
used is assumed to be constant, then the force applied to the ions will be dependant on 
a. This is because the lattice size will increase with increasing a and, therefore, so will 
the spacial area of the beam required to impart a force onto the ions. It is, therefore, 
required to express this force as a function of a. The intensity of a beam required to 
provide a state dependant force, F, is given by [35] 

_ 3FAAJ.,, 

where A is the detuning of the laser from resonance, A is the wavelength of the laser, 
Isat is the saturation intensity of the ion and 7 is 2n times the transition linewidth. The 
power of a laser beam is given by 

P = loa (29) 

where a is the spacial area to which the beam is focused. In this work the beam 
is assumed to be focused to form a light sheet across the array with an area given by 
0- = {ng — I) AW = {ug — l)kA<yW, where n^ is the number of trapping sites (or polygons) 
along one side of the array and W is the width of the light sheet. By using equations 



28] and |29] the force applied to the ions by a laser power, P, can be expressed as 
2TrhP-f'^ 



(30) 



3aAXIsat 
The form of Eq for the case of lasers applying the state dependant force can now 



be found by using equations 08, 30 and the general error equation 27 yielding 



_4:V2 7r^h^ Mkl-f^m''P''a{n+\) 

It follows that both the Kgim and simulation error Eq will decrease with increasing 
laser detuning, A. Therefore, the optimum detuning, for a given laser power and a, 
corresponds to the lowest detuning which provides the required Kgim- The optimum 
detuning to achieve the lowest simulation error for ^^^Yb"^ will be discussed in section 

Magnetic fields can also be used to provide the state dependant force, F, and is 
given by |3] 

'^^d^co/a^^A (32) 
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where u = •ygbi is the position dependant spin resonance frequency with 7^ = e/nie 
the gyromagnetic ratio and i is the x, y or z direction. The magnetic field gradient 
b is assumed to arise from a magnetic field of the form B = Bq + hi, where Bq is a 
constant magnetic field offset. From this, the state dependant force, F--, produced from 
a magnetic field gradient, b~-, is found to be 

„ heb: 



2m, 



(33) 



e 



where rUe and e is the mass and charge of an electron respectively. If one assumes the 
magnetic field is created by a current carrying wire located on the surface of a polygon, 
at a distance a from the centre of the polygon and making an angle of 45° with respect 
to the X-axis then the magnetic field gradient will be of the form 

Here /iq is the permeability of free space, / is the current fiowing through the wire and 
r'^ is the distance squared of the ion from the current carrying wire and is equal to 
r^ + a^ where r is the ion height. We assume, for simplicity, that the distance a scales 
linearly with a with a constant of proportionality of ka in order to keep the angle of r' 
to the x-z plane, 6, independent of a. As the ion height is known to scale linearly with 



a, from equation 15, it is possible to express the magnetic field gradient along r' as 



The component of this magnetic field gradient in the x-z plane can now be shown to be 

_ fipi cos e 

The form of Kgim for the case of magnetic field gradients applying the state 



dependant force can be found by using equations 18, 33, 36 and the general error 



equation [27] yielding 

_ V2 n/XQ fcX/^cos^gM(n+|) 
^-- GAn^mle ,l^{kl + klfa ' ^''^ 



Equations 31 and 37 show that the quantum simulation error is proportional to 
a for a state dependant force created by a laser beam and proportional to a~^ for a 
magnetic field gradient created by current carrying wires. For the case of laser beams 
the a scaling implies that as a is increased (yielding larger K^jm values and geometries 



as shown in section 4.3) the quantum simulation error will rise and, therefore, provide 
an upper limit on the value of a. For the magnetic field gradient case the upper limit on 
a comes from the current creating the gradient. While the a scaling for the quantum 
simulation error using magnetic field gradients implies that a larger a is advantageous, 
the current required to achieve a given magnetic field gradient scales as a^ as can be 
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deduced from equation [35j The maximum current that one can apply to the lattice 
therefore provides an upper limit for a. 

It is also interesting at this point to note the different scaling with a of the laser 



and magnetic field gradient forces given in equations 30 and 32 The laser force can be 
seen to be oc a~^ as it is a function of the inverse of laser sheet cross section, a, which 
is oc a^. The magnetic gradient force, on the other hand, is oc a~^ as it is a function 
of the magnetic field gradient br' which is oc a^^ due to r' having a linear relationship 
with a in the geometry considered. 

6.3. Spontaneous emission 

When applying the state dependant force to the ions using a laser beam additional 
decoherence will occur via spontaneous emission. The spontaneous emission rate is 
given by [36l [37] 




S = ^l— + (3J 

where 7 is 27r times the linewidth in Hz, A is the laser detuning from resonance, Aj^ is 
the fine structure splitting of the ion (fa 100 THz for ^^^Yb"*") and g is the single photon 
Rabi frequency given by 



9 = ^\lwf^- (39) 

Here, Iq is the laser intensity and Igat is the saturation intensity of the ion. It is possible 
to express the single photon Rabi frequency in terms of the laser power, P, by using 



equation 29 giving 



It is now possible to describe an additional parameter, Lgim, which describes the 
ratio of interaction rate to the spontaneous emission rate as 

Lsim = ^ (41) 

where 

Ts = l. (42) 

It has been shown that the detuning which minimises the effect of spontaneous 
emission is ~ 33 THz for ^^^Yb"*" [3H]- It, therefore, follows that the value of L^jm will be 
maximised with this detuning. This additional parameter is analogous to the parameter 
Ksim in equation [3] and is also required to be greater than unity, just like the original 
Ksim, when considering a state dependent force created using laser light. If magnetic 
field gradients are to be used to apply the state dependent force then L^j^ is no longer 
relevant. 
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6.4- Other considerations 

It is important to note here that an increase in a will increase the time taken for ion-ion 
interactions to take place in optimised lattice structures, as we will show in equation 



44 Equation |5] gives an expression for the time taken for an ion-ion interaction to 
occur in any given fixed lattice structure. This can be expressed for optimised lattice 
structures by including the expressions for the ion-ion separation (polygon separation, 



A) and secular frequency, u, from equations 16 and 18 respectively, to yield 



^ 2 k^F^m^a ^ ' 

An expression to give the interaction time in optimised lattices as a function of a 



can now be arrived at by using equations H3^ and 30 



9 eho Vteomns - l)kAwfA'\'ll,a 
^ 8 nh k^^^^P^m? ' ^ ' 

Equation |44] clearly shows that as a is increased the time taken for an ion-ion interaction 
will increase and, so, it may be preferable to limit the magnitude of a after taking into 
account the effects on ¥^sim- A similar equation can also be derived for the use with 
magnetic field gradients. It is also important to note here that increasing the laser power, 
P, will increase the value of Ijsim as the spontaneous emission rate is proportional to 



P whereas the coupling rate is proportional to P^ as shown in equations 38 and 44 
respectively. 

With the use of the equations derived in this section, optimal geometries can be 
calculated given certain experimental parameters and will be described in the following 
section. 

7. Example case study 

In this section, an example case will be presented to show how a 2D lattice for the use 
in quantum simulations can be designed using the work in this paper, whose successful 
operation is within reach of current technology when using lasers. We then go on to 
show that while magnetic field gradients can be used to create a state dependant force 
their use may be more challenging. 



From equation [23] we can see that there is a unique voltage for a given mass, 
m, lattice type and stability parameter, q. For this example case we choose a lattice 
comprised of square type unit cells with 9 trapping sites for ^^^Yb+ ions with a stability 
parameter q = 0.5. Using these parameters the voltage can now be determined by 
calculating k^.. K can be found by plotting the ion height of an optimised lattice against 
a, as shown in figure IS] (a), and finding the gradient of this linear relationship. For this 



example case k^ ~ 98 mV s . Using this result and equation 23 we find the unique 



voltage to be ~ 34 V where rjgeo has been calculated to be ~ 0.145. 
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a 


R [fim] 


A [fim] 


A [THz] 


P [W] 


M 


^sim 


^sim 


Eo 


0.3 


14 


52 


33 


6.5 


1 


35 


1.5 


0.25 



Table 9. Table summarising the parameters for a 3 by 3 square type unit cell lattice 
at cryogenic temperature as shown in the example case study 



The a dependant polygon radius and separation can be calculated by producing 



the corresponding graphs in figure [T0| using the method described in section |4.3[ for 
the lattice type and ion to be used. For this example case the optimum radius and 
separation of the polygons in terms of a is ~ 45a fim and ~ 174q; fim respectively. 

The next step is to choose a laser detuning from resonance and a maximum 
acceptable error, £"0, to be used. For ^^^Yb"*", as explained in section O, the optimum 
detuning is ~ 33 THz giving a wavelength of ^ 355 nm which is therefore used in this 
example case together with a maximum acceptable error of 0.25. Note that in this 
example case the array is assumed to be at cryogenic temperatures which correspond to 
an electric field noise density three orders of magnitude less than at room temperature. 



We can see from equation 31 that a needs to be minimised in order to keep the quantum 
simulation error low. The minimum a is determined by the lowest ion height one can 
easily achieve which for this example case we choose to be equal to 30 fim. For this case 
we find a ~ 0.3. Having determined a we find the optimum radius and separation to 
be equal to ~ 14 /im and ~ 52 fim respectively. To calculate the laser power required 



equation 31 should be set to the maximum acceptable error and solved for the laser 
power, which is found to be ~ 6.5 W assuming the ions to be cooled to rl <^ 1 and 
M = 1 (one side average observable). These conditions provide a coupling rate, J, of 
~ 530 Hz with a /3 ~ 2.8x10"^. This laser power can, for example, be achieved using 
a commercially available diode pumped solid state (Coherent Paladin range) or fibre 
(Coherent Talisker range) laser system. Table M summarises all parameters for this 
example case study. 



Figure 14 shows the effect a change of a and laser power, P, has on the quantum 



simulation error (solid curves), Ksim (dashed curves) and Lsim (dotted lines). The 
cross corresponds to the 2D lattice designed in this example case which represents the 
optimum case in terms of achieving the highest Kgim and Lgim for a maximum quantum 
simulation error of 0.25. We note that the main limitations in achieving lower quantum 
simulation errors stem from the lowest achievable ion height (lowest a) and magnitude 
of electric field noise density. Figure 14 also shows that higher values of Ksim and Lsim 
can be achieved with any given geometry (given by the value of a) by simply increasing 
the power of the laser. However, this can only be achieved at the expense of higher 
quantum simulation errors. 

State dependant forces can also be created using magnetic field gradients as 
described in section 6.2 Figure 15 shows the Ksim (solid curves) and the quantum 



simulation error (dashed curves) as a function of the magnetic field gradient, b, and a 
for traps operated at cryogenic temperature. Here we use ^'''^Yb^ ions in a three by 
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a 



Figure 14. Graph showing the quantum simulation error (solid curves), K^^^ (dashed 
curves) and J^sim (dotted lines) for a three by three square type unit cell lattice with 
^^^Yb+ ions as a function of a and laser power. Here traps are operated at cryogenic 
temperature. The cross indicates the example case. 



three square unit cell array. As described in section 6^, Eq oc a ^ indicating that a 
large a is advantageous. The limit on the maximum a is dependant on the maximum 



current one can apply to the geometry. Using equation |35] it is possible to calculate the 
current required to create a desired magnetic field gradient. In order to illustrate the 
magnitude of currents required we assume ka = kr, which will result in an angle 6 = 45° 
(refer to section 6.2 for more information). We also set a to ~ 0.3, determined by the 



lowest achievable ion height which, for illustration purposes, we have set to 30 jj,m. The 
reason for choosing the minimum a value can be seen when considering equation 35 
which clearly shows that, for a given magnetic field gradient, J oc a^. In the magnetic 
field gradient case, the chosen a sets K^j^ and Eq. For this case, again, we choose M = 
1, n <C 1, Ksim = 2 and Eq ^ 0.01 which requires a magnetic field gradient of ~ 33,000 



Tm~^ and is indicated by the cross on figure 15 This is achievable with a current of 
^ 1,200 A yielding a coupling rate, J, of ^ 240 Hz and a /3 ^ 2.8x10"^. From this 
simple example case one can conclude that using magnetic field gradients to provide 
state dependant forces for the use in quantum simulations using the methods and trap 
designs shown in this work is quite challenging. However, geometries trapping ions in 
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chains allow for sizeable magnetic gradient induced couplings |39] and, so, a detailed 
investigation into optimising the wires used for producing magnetic field gradients in 
the geometries discussed in this work could improve results. 
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Figure 15. Graph showing the quantum simulation error (dashed curves) and the 
Ksim (soHd curves) of a three by three square type unit cell lattice with ^^^Yb+ ions 
as a function of a and magnetic field gradient. Here traps are operated at cryogenic 
temperatures. 



8. Conclusion 



Two-dimensional arrays of surface ion traps have the potential to provide a technology 
with which quantum simulations can be performed. In order for ion traps to be used 
successfully for this purpose a greater understanding is required of how the various 
geometry parameters affect the ions trapped above them. Throughout this work square, 
hexagonal and centre rectangular unit cell arrays of microtraps have been modelled 
in the gap-less plane approximation using the Biot-Savart like law in electrostatics 
|30j . Decoherence due to motional heating [18] was then compared to the ion-ion 
interaction [21] to provide a ratio used to describe how much faster an ion-ion interaction 
occurs in comparison to the motional decoherence, K^jm- This work investigates how 
various parameters in the array can be adjusted in order to optimise the device's 
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ability to perform quantum simulations and shows how the interactions can be made as 
homogeneous as possible over the device. It has been shown how the homogeneity of the 
Ksim across an array can be altered by varying the distance from the outer polygon to 
the edge of the rf electrode. The distance required to maximise the Ksim homogeneity is 
also shown to vary as a function of the total size of the lattice. The number of polygon 
sides, n, required to maximise Kgim has also been found. 

We have shown that the K^j^ of a given lattice geometry can be maximised by 
reducing the value of a. However, as a reduces so does the trap depth. This results in 
the conclusion that the maximum K^j^ of a geometry can be achieved by reducing the 
value of a until the trap depth reaches a reasonable minimum value. 

Using this information, optimal geometries as a function of a are presented. This 
has been achieved by finding the relationships of polygon separation and radius with 
a for optimal geometries. It was found that, for these optimal geometries, Ksim scales 
as a^. The individual polygon separation and radius were found to posses a linear 
relationship with a and, therefore, larger geometries have been found to produce larger 
values of Ksim- Therefore, the optimal lattice geometry is dependent solely on the value 
of a for a given ion mass and number of trapping sites in the array. 

We presented a case study for determining an optimum geometry consisting of 9 
trapping sites arranged into square type unit cells for ^^^Yb"^ ions. We showed how 
the value of a can be chosen (and, therefore, the geometry dimensions) by taking into 
account the laser power (or static magnetic field gradient) required to produce a state 
dependant force acting on the ions, the Ksim and the error on the simulation. From 
this it has been found that to carry out quantum simulations with reasonable Ksim and 
error values it is preferable to use traps held at cryogenic temperatures as this reduces 
decoherence due to heating effects on the ions. Other methods known to significantly 
decrease the anomalous heating rate include pulsed laser electrode cleaning [10] and 
Argon-ion beam electrode cleaning [H]. 

The scaling of anomalous heating has not yet been fully understood and is thought 
to posses a dependence with the ion height, r, of between r~^ and r~^. In this work 
we have used r~^ however, if the anomalous heating is found to posses a relationship 
with ion height nearer r~^ then the equations in this work can be adjusted and this is 



discussed in more detail in [Appendix B 



The relationships between lattice size and a with the polygon radii and separation 
obtained using the method described in this work will allow for the construction of 
two-dimensional surface trap lattice arrays with high ratios of ion-ion interaction rates 
to decoherence rates, providing a system which could be used to perform quantum 
simulations. 
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Appendix A. 



The simulation error in equation 27 can be adjusted to take into account the heating 
of the ions during a simulation. If this occurs the mean radial phonon number n will 
become time dependent, n{t). It, therefore, follows that the error will also become time 
dependent and is given by 

, , F^Minit) + \) , , 

The time dependent mean radial phonon number will be a function of the heating 
rate n, the interaction time, Tj, and the initial mean phonon number, hq and is given 
by 

n{Tj) = no + hTj (A.2) 

yielding an error given by 

It should be noted here that hTj = 1/ Kgim and so the error can be expressed as 

EoiTj) ^ \ ^- '\ (A.4) 



Equation A.4 shows that as Ksim increases the error will tend towards that in 



equation 27 in section 6.2 for the case of n <^ 1. This is because higher values of Kgi^ 
mean that less heating takes place during an interaction until the mean radial phonon 
number can be approximated as constant. 

Appendix B. 

In the work presented we have used a motional heating rate, n, which has an r~^ scaling, 
where r is the ion height. However the scaling of this anomalous heating has not yet 
been fully understood and so it is conceivable that a different scaling will be required 
to describe the effect. For this reason it is the aim of this appendix to outline the steps 
and main expressions required to obtain an optimised two-dimensional ion trap array 
with an anomalous heating rate which posses an arbitrary scaling with the ion height 

For the case of optimising the homogeneity of Kgim across an array the results 



shown in section 4.1 will hold for any scaling of the heating rate with ion height. This is 



Optimisation of two-dimensional ion trap arrays for quantum simulation 29 

because this optimisation aims to give each trapping site in the array the same properties 
and this is independent of the heating rate. The same apphes to the optimisation of the 



number of polygon sides described in section 4.2 



Equation [19] describing Kgim for an optimised geometry can be altered to take into 
account a different scaling of the heating rate with ion height. For an arbitrary scaling 
r~^ this equation can be expressed as 



It can be seen in equation B.l that the Kgim is proportional to a^^ ^'. The values 



of kr, kA and ku (described in section 4.3) are independent of value of x and, so, the 



optimum geometry for a given a will be the same regardless of the scaling of the heating 
rate with ion height. Using this procedure optimum geometries can be computed for an 
arbitrary scaling of the anomalous heating with the ion height, r~^ . 
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